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Sununary 

Logical analysis of a box \ving necessitates the cillowa/nce 
for the contribution of the drag spars to the torsiona.1 strength 
of the structure, 

A rigorous analysis is available in the use of the Method 
of Least Work. 

The best logical method of analysis is that applying 

Prandtl's Hembra.ne Analogy, in the form 

St - ^ 
^ - 2A • 

The results so obtained vary by a negligible amount from those 
obtained by the rigorous method. 

The stresses in the members of a box wing should be ca.1- 
culated by the membrane analogy method, but should be subject 
to verification by test before being used in design. 

Scope 

The scope of this pa-per is the a.nalysis of the elements of 

aw conventional type of box wing under a torsional load. This 

winp: has as a primary structure two ;70oden box or I beams, their 
♦Thesis submitted in partial fulfillment of the requirements for 
the degree of Engineer in L'Iecha«nical Engineering Aeronautics, 
Stanford University, 1930. 
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maximum moments of inertia being about axes parallel, or nearly 
so, to the wind chord. They are connected by a plywood skin, 
forming the wing covering, in such a way that a cross section 
parallel to the plane of symmetry of the airplan-6 forms a 
rectangle or a trapezoid* Figure 1 shov/s a typical box-wing 
cross section, the tv/o spcirs proper being box beams, and the 
skin covering being plywood which forms part of i:he airfoil 
section. None of the formulas covered in this report take ac- 
count of the curvature which in practically all cases is present 
in either the top skin or both top and bottom. It is believed 
that this factor is so unimportant that it would not be v/orth 
while to introduce the complexity attendant upon its consider-- 
at ion into the relatively simple formulas obtained when the 
curvature is neglected* 

It is shown in this paper that the analysis of a box wing 
by rational methods results in the computation of much lower 
stresses in the various members of the box than are obtained 
when present design procedure is followed. As a means of 
simplifying the computations necessary in such an analysis, it 
will also be demonstrated that approximate methods, involving 
some reasonable assumptions, are available. Assuming the 
validity of this statement, which shall be proven subsequently, 
the value of this paper v/ill lie in the application by the 
designer of its methods to reduce the material necessary to 
carry a given wing load; and by so doing, he will decrease the 
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weight of the wing structure. 

Definitions and Nomencla-ture 

Lift spar : A spa.r formed "by two chord members connected hy 
a web member, the chord members lying in a plane a,pproximately 
perpendicula-r to the wing chord. 

Drag spar : A spa.r formed by two chord members connected by 
a web member, the chord members lying in a plane approxima,t ely 
parallel to the wing chord. 

Ela.stic centrum : A point in the wing structure in such a 
position thsit if the line of action of an imposed loo.d pa.sses 
through it, the load will cause no rotation of the cross section. 

Included statically determinate structure : The part of a 
redundant structure which remains when enough of the redundant 
factors have been eliminated to make the remaining structure 
statically determinate. 

Beam force : A force parallel to the intersection of the 
plane of the lift spe-r web and the plane of sjamuetry. 

Chord force : A force pa^rallel to the plane which bisects 
the dihedra.1 angle, or the distance, between the pla^nes of the 
two drag spars or trusses. 

Conventions for Signs^: 

Forces: An upward beain force is positive. A rearward 
chord force is positive. 

Moments: The torsional moment on the wing is a pitching 
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moment; therefore, we will consider a torsional moment positive 
which tends to increase the angle of attack. 

Theory 

The problem covered in this discussion is the division of 
torsional load between the -^^arious elements of a box wing. The 
necessity for nev/ methods of analysis lies in the fact that the 
design rules of the Department of Commerce at present assiime, in 
effect, that the resistance to torsion of the primary structure 
of a box wing is confined to the bending strength of the lift 
spars, that any lift load applied to the wing is divided between 
the lift spars in inverse ratio to the distance of the load from 
them, and that any twist on the wing is carried as a pair of 
equal and opposite beam loads^ This method will result in the 
application of very severe loads upon the rear spar in the' 
required design conditions of low angle of attack and nose dive. 
With most airfoils in common use, the center of pressure of a 
positive air load on the wing moves forward as the angle of 
attack is increased from a position aft of the trailing edge at 
zero or negative angles to a maximum forward position at about 
thirty per cent of the chord aft of the leading edge as the 
attitude corresponding to maximum lift is attained. In a nose 
dive, the forces acting upon the airplane are a dovm load on the 
tail, an up load on the rear spar, and a down load on the front 
spar ^ a condition corresponding to a resultant load acting a 
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chord length or more aft of the trailing edge. If v;e now resolve 
our load into an equal load acting at the elastic centrum and a 
torsion about the elastic centrum, we see that the torsional load 
is exceedingly severe in nose dive and this condition of flight 
Vvfill in a great many cases be the critical one for the reax spar. 
In all except very unconventional designs, the low angle of 
attack condition will be the critical condition of flight for 
the rCcLT spar when nose dive is not; and in this condition as 
well, the resultant load is well aft of the elastic centrum, 
with a consequent high value of the torsional moment. 

Due to the fact that the airfoil section limits the heights 
of the spars, with the rear spar as a general rule being the more 
shallo'.7 of the tv/o, for a given value of the load the strength 
weight ratio of the rear spar tends to be smaller than that of 
the front spar. The saving in weight v/ill be a material gain if 
it can be proven by a logical method that the actual loads in 
the rear spar v^ill be smaller than those calculated by the de- 
sign rules of the Department of Commerce. 

The limitation of space prevents the consideration of more 
than one type of v/ing. For that reason a full cantilever, all- 
wood, structure will be chosen, with the lift spars formed by 
two box beams, and the drag spar web formed by a plywood skin* 
A cross section of a typical wing is shown in Figure 1. Assume 
a torsional moment M applied at the elastic centrum. If the 
strength of the drag spars in torsion is neglected, the moment 
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is resolved into the loads Wj, and w^^ acting upon the lift 
spars, a.nd equal to ^ . Now assume that a certain portion 
of the moment, M^, is carried by the drag spars* Then w^, 
and are obtained from the approximate formula 

And also Mi 

A uniform load of Wj, on the front spar, shown in Fig^are 1, 
will result in a bending moment on the spar, resisted by 
compressive stresses in chord member A and tensile stresses 
in chord member B. The load w^ in the upper drag spar, by 
the sane reasoning^ causes compressive stresses m member C 
and tensile stresses in member A* The lower spar, under the 
load w^j is subjected to compressive stresses in member B 
and tensile stresses in member D. Lastly^ the rear spar^ 
under its load w^^, is subjected to compressive stresses in 
member D and tensile stresses in member C. Member A is 
then subjected to compression from w^, and tension from w^; 
member B to tension from w„ and compression from w, , 
All four chord members of the lift spars are not only with- 
standing a smaller load, but actually are subject in addition 
to stresses of an opposite sign which reduce their net stresse 
still further. 
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Referring to Figure 1 again, it is seen that under a 
torsional load the structure is redundant - we have four 
memhers and but three equations of equilibrium. The redun- 
dancy necessitates, for a rigorous solution, the use of one 
of the methods of consistent deformations, such as the method 
of least work. If rigor is not mandatory, certain simplifying 
assumptions may be applied, such as the neglect of the work 
done in shear, or the assumption of an undistorted cross section 
after the loading. 

Four different basic principles are applied to the problem 
at hand to determine the eight formulas used here. The first 
neglects the strength of the drag spars; upon this, present 
design rules are based. The Theorem of Least Work generates 
the formula of that narne and also the Inverse Ratio Method. 
The simplifying modifications of the Theorem of Least Work are 
responsible for the simplified Method of Least Work, the 
trapezoidal method, Niles' method, and Burgess" Moment of 
Inertia method. The Membrane Analogy is the basis for the 
formula of that name. It is reasonable to expect some corre^ 
spondence between the results of methods developed from the 
same groundwork, and this expectation is realized in the work 
done. 

The formulas to be derived in the remainder of the report 
will be applied to a sample spar, shown in Figure 2, in order 
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to obtain a comparison of the results. The length of this spar 
is 200 in., and the type is a full cantilever; the dimensions 
of the cross section are shown in the figure. 

Present desi.Q:n practice .- The strength of the wing in 
torsion, neglecting the drag spar contribution, is equal to the 
bending strength of the lift spars- A torsional load will then 
be resolved into equal and opposite loads on the lift spars, 
and zero loads on the drag spars. Applying a torsion of unity 
per inch of span on the spar in Figure 2, the following values 
of the running beai'/i loads are obtained^ 

Wji = 0.0286 lb. /in. 

v/„ ^ -0^0S86 Ib./in, 

Jci 

= 0.0 lb. /in. 

w, = 0.0 lb. /in. 
L 

The computations for this, and the remainder of the applications 
of formulae, to the sample spar, v/ill be found in the appendix. 

The Method of Least Work.- The method of least v/ork was 
developed and proven rigorously by Castigliano; it states that 
the internal work done in a redundant structure by the appli- 
cation of external loads will be the least possible consistent 
with equilibrium- The deriva.tion and proof of this theorem 
will not be given here, since any textbook on elasticity or 
*See Andrew^s "Elastic Stresses in Stractures." 
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advanced mechanics of materials contains the development in full. 
The application of this method to the problem at hand will take 
the following form: an expression for the total v^ork done in 
the structure, in the terms of one of the unknown loads, will be 
set up, differentiated with respect to the unknown, set equal to 
zero, and solved. Since there are four unknowns and three 
equations of equilibrium, all exc'ept one of the unknowns may be 
eliminated from the expression for work, and one differentiation 
will serve to determine the solution* 

The bending work done on a beam of constant cross section 
is easily expressed as 



M = bending moment of external loads 

L = length of beam 

E = modulus of elasticity 

I = moment of inertia of cross section 

X = distance along span. 



This does not express the total work done in the structure, 
however, since the shearing work done in the beam must also be 
considered. The following development for the internal work 
done in shear on a box beam is taken from N.A.C.A, Technical 
Report No. 180, Reflection of Bea^s with Special Reference to 




^ Mfdx* 
2EI 



where 



the internal work done in bending 



*Spofford^s "Theory of Structures." 
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Shear Deformations," b3'- J. A. Newlin and G. V/. Trailer. On page 
18 of the report, vie find 



- V 



t ^ydy + t jdy 



aw^ =gtdyeu^=t,dy<u2w^ (3) 

+ 2 / ^ ^ t dydx (4) 

o y 2F ^ 

where L = seniispan, or len^tli of bean (in inches) 

F = shearing modulus (for spruce, l/l5 E) 

I = monent of inertia of cross section 

K = distance fror.i neutral axis to flange 
1 

K = distance from neutral axis to extrene fiber 

t = web thickness (inches) 
1 



2 

q 



t^ = flange width (inches) 

- unit sh earring stress (ib./sq.in.) 

■ unit shearing stress in flange (ib./sq.in. ) 

= unit shearing stress in web (ib./sq.in-) 

V = total shear (lb.) 

d = 2 ( syrometricaJ beam) 

y = distance from neutral axis 

X = distance a.long span 
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It -V 2^ ^ 



(5) 



It 



t^ydy + /^i t^ydy 

\ y 



It 



O 2 1 <5 1 



(6) 



From 4 



4 



V^t 



8FI' 



2F ~ 2F 4 



t_ , K _ -trS 



8 _Jr /V 2 
1 41^ 



15 



(7) 



Let 



a 



8FI' 



8K 



15 



(8) 



and 



•Wc, = / a dx 



2"S 



o 



(9) 



o 2F i 



t K, 

2F o 



Pt^2 



dy. O-O) 



Eq. 10 = 



2FI^t. 



2 1 



(11) 



13 



Let 



7 = 
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t ..""K 



■t t K ^ o ^ = 



Then 



W_ = 2(a + 7) V^djc 
o o 



(12) 
(13) 
(14) 



For a cantilever boani, uniformly loaded with w lb. /in., 

V = wx ; (15) 



W„ ^ 3(a + 7) / w^x^dx = | (a + 7) w^L' 
o n 3 



(16) 



Equations for the internal work in a v;ing whose spars are 
not uniform in cross section are seldom expressible as func- 
tions of X. In such cases, a unit length of span will toe 
treated as a uniforn section, the value of the internal work 
on that unit length found, and the resulting equation treated 
exactly as the one here obtained to get the load division at 
the one point;, the process is repeated until a curve of load 
division is defined, and then the running loe.ds on each com- 
ponent of the box are known* The values of a and 7 as 
determined by equations (s) and (12), are for synnietrical beans 
only. For beams having unequal chord members, only a very small 
error is introduced by using a and 7 determined for an 
equivalent symmetrica.l beam having the same over-all height, 
width, and the Scime moment of inertia. This approximation will 
be necessary in almost all cases for computations concerning 
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the drag spars. 

Having the expression for the bending work done, and the 
shear work from equation (16)^ in terns of the running load on 
the component of the box being considered, the work done will 
be expressed in terms of one variable by substitution from the 
equations of equilibrium, obtained by inspection from Figure !• 

Wy cos ri - Wj^ = 0 (17) 

Mom. at B 

Wy h^ cos T) + dWp^ = 1.1 =1 (18) 
(for moment of unity). 
Mom, at D 

h^ cos r| 4- dWj, = 1 (19) 

From these three equations, all except one of the unknowns in 
the total work equation for the spar may be eliminated; having^ 
the total work, then in the form 

= f(w^) (20) 
differentiate W^j^ with respect to w^; 

^ = f Cwx) 0. (21) 

This serves to determine w^, and from the equations of 
equilibrium the remaining three unknowns may be found. 

Applying this method to the spar of Fignire 3, the fol- 
lowing results are obtained by the calculations shown in the 
appendix: 
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= 0.0168 lb, /in. 




= -0.0134 lb. /in. 




= 0.0590 lb. /in. 




^- -0.0589 lb. /in. 



An examination of the work in the appendix under this 
method illustrates the reason for not using the theorem of 
least work more often. The calculations are involved and 
tedious, with a great ma.ny chances for error, and no inde- 
pendent check. It does serve the useful purpose, in a paper 
of this kind, of constituting a check with which the other 
approximate methods may be compared. 

Simplified Method of Least Work .- The example used for 
the computation by the least work method showed that the per- 
centage of the total v/ork done in shear was small. As a means 
of obtaining a simpler solution, the shea.ring Y^ork will be 
neglected as an approximation, and the load division will be 
found on the assumption that all of the work done is done in 
bending. The method for this calculation is exactly the same 
as for the complete least work method, the only difference 
being that the v\^ork done in the member has only one term in it 
instead of two. The expression for totaJ v/ork is again ob- 
tained in terms of one variable by the sej.ie equa::ions of equi- 
librium; the equation is differentiated with respect to the 
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variable, equa,ted to zero, and solved. The results on the 
sample spar are: 





0.0144 Ib./in* 




-0*0103 To. /in. 


= 


0.0712 Ibi/in. 




-0.0710 lb. /in. 



The difference between these results ajid those obto.ined 
by the rigorous method is due to the fact tha,t the shearing 
work done in the different elenents of the box is not a con- 
stant percentage of the tota.l work in eo.ch elenent. In the 
front spar, the percentage of shea-ring work is in the 

rear spar, 4.4^; in the top a.nd bottom spars, 38.1^. This 
is not unexpected, since the proportions of the elements are 
80 dissimilar. As the spar becomes deeper a,nd thinner, the 
bending work under a given load becomes less, while the 
shearing work is more (for a given cross-sectiona.l area). 
The variation in the running loads obtained from those ob- 
tained by the rigorous solution is, for the lift spars, -14.3^ 
in the front and -23.lfo in the rear; for the drag spa.rs, 20* 7^. 

Solution of a tra.pezoid >>- Referring to Fi^^ure 1, it is 
seen as before that the three equa-tions of equilibrium a.re 
insufficient to determine the four unknom loads. To elimina.te 
the redundancy, the assumption shall be ma-dc that the cross 
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section suffers no distortion after a torsional loading; and 
as a simplifying assumption, consider the deflections of the 
various elements of the box as inversely proportional to their 
moments of inertia, which corresponds to the previous simpli- 
fying assumption of negligible shearing work. 

If the cross section suffers no distortion after loading, 
the change in slope of two sides of the box ma.y be equated; 
and the change in slope will be expressed by the equation 

An - 6p - - 5tj sin r\ ^ 6^ cos r\ ^ 

where 

0 = torsiona.1 angle of twist 
6 = beam deflection 

r\ = angle between sloping dr a.g spar and 
spar opposite 

6 = (2) 

where K is determined by the elastic curve of the beam, and 

will be the same for all elements of the box, on the assumption 

that the distribution of the drag load is similar to the lift 

load distribution. By substituting in (l), we find 
]_ r K Hp K Wpj Z v/rj sin r\ 



r 



^ 1 

ii 



K COST)' K 



(3) 
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This equa.tion, plus the three equations of equilibrium given 
before, will be sufficient to determine the values of the dif*- 
ferent loads. The solution of the set of four simultaneous 
equations in a general form is very cumbersome and should not 
be attemptedc The process is greatly simplified by substituting 
numerical values in the simultaneous equations and then solving 
by any one of the standard algebraic methods. Such a solution, 
applying the constants of the sample spar, results in 

Wj, = 0.0140 lb. /in. 

Wj^ = -0.0098 lb. /in. 

Wy - 0.0731 lb. /in. 

= -0.0729 lb. /in. 

The accuracy of this method of solution is poorer than the 
simplified least work results, the differences in the lift spar 
loads being -16.7fc» in the front and --26,9fo in the rear; the 
drag spar difference being 23.8^. 

Uiles' Method of Load Division .- The basic assumption under- 
lying a method developed by A. S. Niles, is that the cross sec- 
tion suffers no distortion during the application of a torsional 
load. As simplifying weapons, we a-lso assume that the drag 
spars are mutually parallel, and are perpendicular to the lift 
spars, and that the shearing work may be neglected. The last 
two of these assumptions simplify the equation equating slope 
increment, to the form: 
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(1) 



d h 

The term h must be approximated here. It seems reasonable to 
use that value of it which is obta.ined at the elastic centrum, 



assumine- a linear variation from h to h. 



This value is 



h = h^ - I (\ - (2) 
A f-urther simplification lies in the fact that Wp and 



as 



are eo.ual, as well as Wy and w^. 



and 



h ' 



They can be expressed 

(3) 
(4) 



where M^^ is the portion of the tots^l moment which is resisted 
by the lift spars, and 1% that resisted by the drag spars. 
Substituting these values in (l), 

r 



1 

d 



K II. 



n 



K 1 



dl, 



d 



1 

h 



K 1% 



K li 



hi 



+ 



a 



hli 



(5) 



= d 



I + 1 



(6) 



This does not give \TOrka.ble results- The fact that the 
section is a trapezoid means that there will be a lift component 
of the load in the slanting drag spar. Therefore, to obtain 
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2V = 0, we must correct the values of w obtained from equa- 
tions (3) and (4) in the follov;ing nianner: 

w.. = _JLi_ 8) 

U COS Tj 

Wp - Wj^ = w«y sin T) (9) 
(d - x) = Mj^ (10) 

These values of wj^, wt^, and w^, the true values, will 
be proportioned so that the total moment on the section is the 
same as that for the first com.putation. 

By substituting the constants of the sample spar, we obtain 
the follov/ing results for unit torsion per inch of span. 

Wj, = 0.0142 lb. /in. 

= -0.0101 lb. /in. 

w^ = 0.0719 lb. /in. 

w = -0.0717 lb. /in. 
L 

These results are relatix^ely close to those obtained from 
the trapezoidal derivation. This is a reasonable indication, 
then, that the assumption of parallel spars, as far as the mo- 
ment division is concerned, involves no major additional error* 
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Burgess ^ Hethod of Load Division .- C. P. Burgess has 
developed* a formula for load division betv/een the various 
elements of a box spar which depends upon the basic assump- 
tion of negligible shearing work in the spars; however, hj 
the same implicit assumption, a,s in Niles* method the results 
are in error by the lift component of loa.d in the slanting drag 
spar, and must be corrected for that. 

The development is as follows: Let s^ be the distance 
of the member x from the elastic centrum. Then the torsional 
rigidity of the member x^, on its resistance to torsional load, 
is I„s^,. The moment of its torsiona.l resista.nce is I^s^^. 

X ^ XX 

Therefore, the load w^ in the member x is expressed by the 
equation 

w — T.-r . .^v- ^ 

% - ^-^zT^ 

The application of this fornula, and the corrections, to 
the Scimple spar gives the following results: 

= 0.0142 lb. /in. 

Wj^ - -0.0101 lb. /in. 

= 0.0719 lb. /in. 

= -0.0717 lb. /in. 

The close agreement between these results and those ob- 
tained through the application of Niles' method would lead one 

to believe that the formulas are simile.r. This is true, a/nd 
♦N.A.C.A. Technical Report ITo. 329, "The Torsional Strength of 
Wings," by C. P. Burgess. 
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the expansion of Burgess' formula into Niles' is given in the 
appendix. For the case in hand, that of a two spar box, either 
formula is equally convenient; Burgess* has the advantage, how- 
ever, of being more easily applied if the structure has more 
than two spars • 

The Membrajie Analogy Hethod*- If we may assume that the 
walls of the box spar are thin with respect to their height, 
we have at our command a fonnula developed from L, Prandtl^s 
membrane analogT«* The derivation of this formula will not be 
given here, but is fully explained in the reference given be- 
low. In its basic form, the equation is 

^ ~ 2At 

where S = shearing stress, 

Q = torque on section, 

A = area enclosed by the centerlines of the sides, 
t = thickness of the side being considered. 
A modification of this formula will be more useful than 
this basic form. If we multiply both sides of the equation 
t, we obtain 

^ - 2A 

where S* is the shear per inch of perineter of the cross 

section. The value of need then be multiplied only by the 

width of the side to obtain the running: load v/. ApDlving this 
•••See Timoshenko and Lesells "Applied Elasticity^* pp. 45 et seq. 
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formula,, we obtain 



and 



w. 



S' = 



F - 



0.001784 



0.01606 



lb. /in. of perimeter 
lb. / in. 



R ~ 



0.01249 



lb. /in. 



u - 



0.0627 



lb. / in. 



Wt = -0.0625 



lb. /in. 



The only error which enters into this calculation lies in 
the variation of the front ejid rear spars from a thin walled 
section. The accuracy, referring to the least work calculation 
again, is excellent; the errors are -4.2^5 in the front and 
-7.2fa in the rear spar, c.\nd 6.3fo in the drag spars. 

Burgess ' Inverse Ratio Ilethod for Load Division . - G . P . 
Burgess has also developed a method called the inverse ratio 
method for determining load, division between the various parts 
of a redundant structure,,* The basic theorem is that the por- 
tion of the imposed load carried by each of the included stat- 
ically determinate structures is inversely proportional to the 
internal work done in the included structure when carrying the 
whole load alone. This theorem is true only when no part of 
one included structure reacts upon any part of any other in- 
cluded structure; when no member of the stricture is common to 
two or more of the included statically determinate systems into 

which the structure is divided; and another error which may 

♦See Airship Design, by C. P. Burgess. 
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become involved in the results of this formula is generated 
when the a,pplied loo.d is not cone entrant od at a single point 
common to the determinate systems but is distributed among 
two or more comon points. For these rea^sons this method is 
not rigorous in all cases. Applying it, we obtain 



Wji = 


0.0166 lb, /in. 




-0.0132 lb. /in. 




0.0595 lb ../in. 




-0.0594 lb. /in. 



The accuraxy of this method, comparing it with the results 
from the theorem of least \7ork, is very good. As a rigorous 
method, it faJls down becajuse the various members interact under 
load. However, in a simpler case than the one at hand, there 
is reason to believe that the a.greenent will become absolute, 
and such a case, one in which the center lines of the spars form 
a recta.ngle, is analysed in the appendix^ The results obtained 
from inverse ratio and least v;ork agree exactly. The reason for 
the a-greement lies in the fa.ct that there is no component of load 
from the dra^g spars entering into the lift spars. Since in 
practice the rea.r spa-r is almost alv/a.ys more shallow than the 
front, this case is a-pparently of aca.demic interest aJone. 

The variation of the results obtained by inverse ratio, 
from those of the least work method, is smaller than for any 
other simplified method. The calcuLations necessitated for 
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this method are more tedious thaji those for the rigorous method 
and for this reason the inverse ratio method is at a disa.dvantage 
in comparison with the membrane analogy method. 

S'unimary of Results 



Method 


lb. / in. 


i 

er- 
ror 


llD./in. 


fo j Wy 

|lt./ia 

ror 


er- 
ror 


lb . / i n . 


P 
er- 
ror 


Present design 


0.0286 


+ 70.2 


-0.0286 


+113.410.0000 


-100.0 


0.0000 


-100.0 


Least \7ork 


0.0168 


0.0 


-0.0154 


0.0 iO. 0590 


0.0 


-0.0589 


0.0 


Simplified 


0.0144 


-14.3 


-0.0103 


-23.1 0.0712 


+ 20.7 


-0.0710 


^0.5 


least v/ork 
















Trapezoidal 


0.0140 


-16.7 


-0.0098 


-26.9 !0. 0731 


+ 23.8 


-0.0729 


+ 23.6 


Niles^ method 


0.0142 


-15 .5 


-0.0101 


-24. 6! 0.0719 


+ 21.9 


-0.0717 


+ 21.7 


Burgess^ mom. 
of I. 

Membrane analogy 


0.0142 


-15.5 


-0.0101 


-24.6:0.0719 

i 


+21.9 


-0.0717 


+ 21.7 


0.0161 


-4.2 


-0.0125 


-7.210.0627 


+6.3 


-0,0625 


+6.1 


Inverse ratio 


0.0165 


-1.2 


-0.0132 


-1.5:0. 0595 

• . . , , u . ■ 


+0.8 


-0.0594 


+0.8 



-~ 1 . .,1 i i : i i I , , 

Values in ta."ble are "beajii loads per inch of span for a torsional moment of- 
unity per inch of span on the spar of Figure 2, 



Discussion of Results 

Present design practice, as the values of the loads in the 
summeory show, is extremely conservative in computing the loads 
in the spars of a, box 'vTing, and at the same time does not pro- 
vide any mea.ns of computing the considere.ble shearing stress in 
the plywood skin. Due to the action of the drag spars, a large 
portion of the torsional load is removed from the lift spa.rs, 
and stresses of a sense opposing the stresses alrea^dy present in 
the chord members of the box are set up. It must be realized, 
at this point, that the ca.lcula-ted stresses detailed in this 
paper are all obtained upon the a.ssumption of a perfect structure, 
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one in which there is no give in the joints, and more particu- 
larly, one in which all the web nenbers, including the skin, 
have no tendency to buckle under a shear load. In practice, 
the skin of the wing will not carry shear without a certain 
tendency to buckle, a.nd for this reason will not maintain the 
theoretical transference of load from web to flange. Under 
these conditions, the a^ctual stresses in the wing will differ 
from those obtained in the rigorous analysis, and will tend 
toward the present design condition. For this reason it is 
not possible to recommend the a^doption of any logical analysis 
without the support of test data. A superficia^l consideration 
of the problem will suffice to show that the internal work in 
the spar under 8. torsional load will be the least when the net 
axial load induced by bending moment in the chord members ap- 
proaches zero^ This condition is that demonstrated in the 
method of least work. 

The application of the method of least work is much too 
cumbersome a means of attack to use when, at the expense of a 
relatively slight loss of accuracy, much simpler methods are 
available. It is probably the simplest rigorous solution 
which can be applied, and for tha.t reason has been used as a 
check against the approximate answers obtained. 

The method of simplified least work has little advantage 
over the rigorous method, since the loss in accuracy is so 
large. It is undoubtedly trae that the labor of computation 
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has been greatly reduced by the omission of the terns of 
shearing; but another approxina^te method is available which 
combines greater accuracy of results with much simpler ca.lcu- 
lations than may be obtained by any simplification of the method 
of least work. For this reason, while the va.lidity of the 
method is recognized^ it has little value. 

The solution of a trapezoidal form of box involves tv/o 
assumptions: first;, that the shearing work done is negligible, 
and second, that the cross section of the wing suffers no dis- 
tortion during rotation. The method of simplified least work 
has a.lready demonstrated the errors attendo.nt upon the first 
assumption. The second has been verified within the limits of 
experimental error on a few occa^sions; its use is, however, 
definitely an a.pproximation. The results of this method demon- 
strate the fact tha.t an additional error does enter into the 
equation when the second assumption, previously mentioned, is 
used. The ma^gnitudes of the inaccuracies obtained by this 
formula, are such as to reduce the value of the method to a 
very small qua.ntity. 

A. S. Niles* development for the treatment of this problem 
may be applied to a trapezoid only if it is assumed that the 
drag spars are parallel to each other. The method is essen- 
tially nothing more than a special case of the trapezoidal solu- 
tion obtained when the lift spar heights are equal, and the 
a.ngle r\ is zero. To apply it to a trapezoid is obviously 
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unreasonable, on that ba.sis. However, if such an application 
is nade, the lift spar loads nust be adjusted to make ZV and 
ZM equal to zero, involving an error from the rational method 
of load division. Including then, as it does, approximations 
in addition to the ones utilized in the previous method, it 
seems logical to expect larger inaccuracies in the results ob- 
tained. This does not occur; the magnitude of the error is 
less for this method than for the trapezoidal solution* Such 
results may not be expected in all cases however, for the rea- 
sons stated above, and therefore this equation should be ranked 
below the tra-pezoidal one in point of accuracy; in regards to 
utility, it is slightly superior, since it does not involve the 
solution of any complicated simultaneous equations. 

C. P. Burgess^ method has as its main asset the ease with 
which it may be applied to a structure composed of more than 
two spars. For the analysis of a two spar wing, it lies on a 
par with Niles^ method, since the tv/o equations are identical* 
The two equations, as far as discussion and results axe con- 
cerned, may be cla^ssed as one. 

The membrane analogy method is by far the best approximate 
method available at the present time for this analysis. It is 
the simplest of all seven of the formulas covered, and the er- 
rors in it are small enough so that it may be used directly as 
a method of computation. The only reason for any error entering 
into the calcula.tions is that the chord members as such do not 
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constitute part of "a thin walled structure;" but the inac- 
curacy clue to this approxination is snail enough to be rela- 
tively unimportant. 

Burgess* inverse ratio nethod, in this case, gives ex-- 
trenely good results. The generality of such an occurrence, 
however, seens doubtful; and it is evident that this nethod 
involves nore labor of conputations than the nethod of least 
work. For this reason, and the fact that there is an approxi- 
nation in the inverse ratio theory, there would be no reason 
for not using least work if sufficient tine and labor were 
available for either nethod. At best, the inverse ratio nethod 
is inferior to the least work nethod in both sinplicity and ac- 
curacy. 
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Appendix 

Sample spar #1 chaxacteristics (see Fig. 2). 
Front spar: 

i(9.0' - 5.5^) + i X i X 9^ = 93.7 + 7.6 = 101.3 in.* 
Rear spar: 

|(7.0^ ~ 4.5^) + i X i X 7^ = 42.0 + 3.6 = 45.6 in.* 

Top and bottom spars: 

Neutral axis = ^^'^ ^ = 13.48 in. from front spar. 

3.50 X 13.48^ + 2.19 X (35.0 - 13.48)^ = 
= 636 + 1014 = 1650 in.* 
Position of elastic centrum: 

35 X 45. 6 _ -IP, 07 
45.6 + 101.3 

(d - x) = 35.00 - 10.87 = 24.13 

\ = y = 4.19 in. 

n = tan"^ A = tan"^ 0.0572 = 3° 17' 

sin Ti = 0.0572; cos r] = 0.998 
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Loads in sample wing, assuming zero drag loads. 
Let 

M = unity per inch of span 
Then Wj, = _wj^ = ^ 

and Wj, = 0.0286 lb. /in. 



= -0.0386 lb. /in, 
Wy = 0.0 lb. /in. 



= 0.0 lb. /in. 



Loads in sample wing - the method of least work. 



From equations (8) and (12), pages 11 and 12, 



^ 8K s 



a = 



and 7 = 



8F 12 
1 

2F I^t^ 



2 



15 



then W. = 2(a + 7) x^dx = -|(a + 7) w^L^ . 

o o o 

From equation on page 9, 

W, = 1^ = ^ (since M = Iwx^for cantilever) 
^ o 2EI 40il 2 



To determine the internal work W total, 

b S 
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Front spar: 



a 



F - 


700 1h /c^.Q in 


= 


2, 75 Ib./sq. in. 


= 


4. 50 Ib./sq. in. 




0.25 Ib./sq. in. 




2.00 Ib./sq. in. 


I = 


101.3 in/ 



a = 



I - 2.75 (4.5* - 

8 X 86,700 X 101.3^ [ 

- I X 4.5^ X 2.75^ 3.75*)| 

- r 1 ^ 

2.81 X 10 986 - 2.75 (411 _ 102 + lO)!- = 

L ^ 

= 2.81 X 108 X 10~^ = 3.04 X 10~® 



2 X 86,700 X 101.3^ X 0.25 L 4 



J 2lx2^ (^^^2 _ 2.75^) 



+ 2 X 0.25 X 3.75^ . ^^^^z _ 2,,^^^) + 

+ ^ X 0.25^ X 2.75^ 1^ = 2.25 X lO"" (2.75 

J 

[20.3 - 7.6]^ + 3.47 [20.3 - 7.6] + 1.3). 



7 = 1.01 X 10 ^ 



Wg = (1.101 + 0.030) X 10 ^x I v.^^ X 200=" = 6.0 w^,^ 
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Rear spar:: 

F = 86,700 Ib./sq.in. 
K ^ = 2.25 Ib./sq. in. 
= 3.50 lb. /sq. in.. 
= 0.25 Ib./sq.in. 
= 1.75 Ib./sq.in. 
I = 45.6 in. * 

a = hl^ /5 X 3.50^ 2 25 r3 5* X 3 5^ X 

8 X 86,700 X 45.5= 2,25L3.5 -3X3.5 X 

X 2.25® + i X 2.25*}| 

a = 12.13 X 10~^° |280 - 2.25 [150 - 41 + 5]| = 

= 12.13 X 23 X 10"^° = 2.79 X 10~^ 

7 = i — {^'^5® X 2.25 (3^53 ^ 2^353^-^ 

2 X 86,700 X 45.6^ X 0.25 4 

^ 1.75 X 0.25 X 2.25' (3^s" ^ 2.25^) ^ 



+ ^ X 0.25^ X 2.25^1 = 11.1 X 10 ^ 

(l72 [12.3 - 5.1]^ + 1.66 Cl2.3 - 5.l] + 0.5) 



y = 1.120 X 10"^ 



W = (1.120 + 0.028) X I X Wj^2 X 200^ = 6.1 w^' 



5 



2 X 200' 

^ 40 X 1.3 X 10^ X 45.6 ^^^'^ 
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Drag spar: 

I = 1652 in.* 
= 0.0625 in. 
tg = 1(1.25 + 1,75) = 1.50 in. 

h = 36.9 in. 
Then if A = area of one flange, 
I = 1652 = 2(3^f A; 
A = 1652 ^ 681 = 2.43 in. 
= 18.45 - ~— = 16.84 in. 
Kg = 18.45 in. 



a = 



1.50 



X 18.45^ - 16.84 [18.45* - 



8 X 86,700 X 1552^ 

- |l8.45^ X 16.84^ + 4 X 16.54*]) 

3 5 J 

a = 7.92 X 10~^^ jl.l4 x 10^ - 16.54 [116000 - 64500 + 

+ 15100] } 

_13 ' - -10 

a = 7.92 X 10 X 0,001 x 10 = 7.92 x 10 - a 

is negligible. 



2 X 86,700 X 1652^ x 0.0625 * 

- 16.84=)^ "^rx^ (18.45^ - 16.84^) + 

+ ^ X 0.0625^ x 15.84^] = 3.34 x lo"^^ 
lo 

[3080 + 8500 + 710] = 4.31 X lO"' 
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: 2.29 w^^ 





W 


= 4.31 X 10"'' X 1 X w^^ X 
Wj^^ X 200^ 




b 


40 X 1.3 X 10^ X 1652 




By 


inspection, from Figure 1 


(1) 




cos T\ - = 0 






cos T) + dYTp = K = 1 


(3) 


1 


cos ri v;., + dv/n = 1 




In 


terms of w 

u 


(4) 


w 


= vjyr COS r\ = 0.998 w,. 

u u 


(5) 




= i(i - COS n wy) = i 






= (0.0286 - ( 


(6) 




- i(l ^. cos ^ " i 



= (0.0286 - 0.257 w^) 

Substituting and collecting, 

Total = X [2.29 + 3.72] (l + 0.998*) +(135.0 + 

+ 6.1) (0.0286 - 0.257 Wy)^+ 

+ (60.8 + 6.0) (0.0286 - 0.1996 w^) 

3Wn 



^ = w,, [1.996 X 2 X 6.01 + 141.1 x 2 x 0.257* + 66.8 X 



X 2 X .1996*] - [66.8 x 2 x 0.1996 X 0.0285 



- 141.1 X 2 X 0.257 X 0.0385] = 0 



47.9 w - 2.88 = 0: 

U 

= 0.0590 lb. /in. 
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= -0.998 X 0.0590 = -0.0589 lb. /in* 

= 0i0286 - 0.1996 X 0.0590 = 0.0168 lb. /in. 

= -0.0286 + 0.257 x 0.0590 = ^0.0134 lb. /in. 

Loads in sample wing - by least workj neglecting sheau. 

By omitting shear terms in total work equation, 

W = (3.72) (1.995) + 50.8 (0.0286 - 0.1996 w..) ^ + 



37.56 - 2.69 = 0; 

= 0.0712 lb. /in. 

Wj^ = 0.0286 - 0.257 X 0-.0712 = -0.0103 Ib./in^ 
= 0.0286 - 0.1996 X 0.0712 = 0.0144 lb-. /in. 
= 0.998 X 0.0712 = -0.0710 lb. /in. 

Loads in sample wing - trapezoidal method. 

By the principles of equilibrium, from Figure 1, 

w cos T) - w =0 
U L 

w h cos ri + w^pd = 1 
U ' ^ 

w„ cos T] + w d = 1 



+ 135.0 (0.0286 - 0.257 w. )^ 




= Wy (2 X 1.996 X 3.72 + 60.8 X 2 x 0.1996^ + 135.0 
X 2 X 0.257^) - (2 X 60.8 X 0-.0286 X 
X 0.1996 - 2 X 135.0 X 0.0286 X 0.257) 
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And by equating angles of deflection, from page 16, (3), 

^ lu ^ -^U -^L 

Substituting numerical values, collecting, 
(iXs) = 0.998 Wy 

(2KS) Wj, = (0.0286 - 0.1996 w^) 
(3X7) vj^ = (0.0286 - 0.257 w^) 

(4>(8) 0.0286 Co. 00987 + 0.02190 _ 0.00003 w^] - 

- 0.111 Co. 0006a w + 0.000603 w ] = 

U L 

(8) K8^ w + 2.21 w_ _ 0.003 w„ _ 0.238 w^^ _ 0.237 w = 0 

F xi U U L 

Substitute 5 , 6 , and 7 in 8' 

(9) (0.0286 - 0.1996 Wy) + 2.21 (0.0286 - 0.257 w^) _ 

- 0.241 - 0.235 = 0 

(9') Wy (-0.1996 - 0.568 - 0.241 - 0.236) = -0.0286 - 0.0 

1,245 = 0.0912; = 0.0731 lb. /in. 

(id) VTj^ = 0.998 X 0.0731 = -0.0729 lb. /in. 

(11) = 0.0286 - 0.0731 x 0.1995 = 0.0140 lb. /in. 

(12) vr = 0.0286 - 0.0731 x 0.257 = -0.0098 lb. /in. 

^ ^ R 

Load division, by JTiles' method. 
1.1" 
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h = 9 _ 2 X 101 3'^+^45 6 = 9 • 0 - 0. 63 = 8. 37 in. 

r?^ ^ - i 2 X 0.000605 \ ^ p^nr. 

- 8. 37 2 b. 00987 + 0.0219^"" ^•'^^^ 

(3) LI,, + lv% = 1.0 

(3)^(2)=(4) = = 0.400; % = -^-^ = 0.600 

(5) w, = 0,600 ^ 0.0717 lb. /in. 

(6) Wtj = 0.0717 0.998 = 0.0719 lb. /in. 

(7) v/p - v/j^ = sin Ti = 0.0572 x 0.0719 = 0.00412 

(8) X + Wj^ (d - x) = 0.400 = 10.87 Wj, + 24.13 vj^ 
(7;->(8)=(9) 10.87 v/p + 24.13 (wp - 0.00412) = 0.400 

(10) 55 -.7p = 0.400 + 0.099 = 0.499; = 0.0142 lb. /in. 

(11) v;^ = v/p - 0.00412 = 0.0142 - 0.00412 = 0.0101 lb. /in. 



Load division by 


Burgess' 


moment 


of inertia 


method* 


Member 


I 


8 


Is 




w 


F.S. 


101.3 


10.9 


1102 


12000 




R.S. 


45.6 


24.1 


1100 


26500 




U.S. 


1652.0 


4.19 


6930 


29050 




L.S. 


1652.0 


4.19 


6930 


29050 
95500 


0.0717 



(1) = 0.0717 -r 0.998 = 0.0719 lb. /in. 

(2) Wj, - Vvj^ = 0.0572 X 0.0719 = 0.00412 



N.A.C.A. Technical Note No. 366 



38 



(3) 10.9 w« + 24.1 Wr, = 1 - 4.19 (0.0719 + 0.0717) = 1 - 



(4) 10.9 Wj, + 24.1 (wj, - 0.00412) = 0.400 

(5) 35 Wj, = 0.400 + 0.099 = 0.499 

(6) Wp, = 0.0142 lb. /in. 

(7) = 0.0142 ~ 0.0041 = 0.0101 lb. /in. 

Proof of identity of Burgess' and Niles' Methods 



~ 0.600 = 0.400 



(2) By Mies 



(1) By Burgess 




(3)By statics, 



% = "^^^ + ^^(d - x) 

% = ^u^^ - y) + \y 



and 



(6) from (4) and (5), ^ 



(5) from (1), 



(4)frora (1), 



Ipx^ + lR(d - x)^ 

Iu(h y) ^ + I^y^ 
Ij,X^ + Ij^(d - x)^ 



(for M = 1) 



Iy(h - y)% I^y^ 



(7) by definition, 



X = 




y = 
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iVL 



(8) (7K6) :^ = 



Ir, + I 



RJ 



lyh 



-.2 



I 

u 



(9) 



h 



^F Ir 
^F + Ir 



1~ 



(10) 



M- 

IT 



•h 



[ If Ir 



h @ 



Q.S.D. 



Loa.d division, by membrane analogy 



S» 

A 



^F 



'R 



Wt = 



'U 



ZA 

35 X i(7 + 9) = 280 sq. in. 

2 ^ 280 ~ 0.001784 lb. /in. of perimeter 
0.001784 X 9.0 = 0.01606 lb. /in. 
0.001784 X 7.0 - 0.01249 lb. /in. 
0.001784 X 35.0 = 0.0625 lb. /in. 
V 0.998 = 0.0627 lb. /in. 
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Load division by inverse ratio. 
From pp.30 et seq. , 

(1) Qj, = 66.7 

(2) Qp^ = 141.0 
^Z) = 6.01 

(4) = S.Ol 

(5) Bp = 35 X i4;l°+-ss.7 = 11.24 in. 
{6) 8^ = 35 - 11.64 = 23.76 in. 



(7) Sy = 0.50 X 8.37 = 4.19 in. 

(8) 8^ = 0.50 X 8.37 = 4.19 in. 

1 Q,-p Q,^ Qh^ 

(9) = 1 + + + = 1 + + 11.1 + 11-1 = 23.7 

0.0422 



23,7 

(10) = 0.0422 X = 0.0199 

(11) qy = 0.0422 X = 0.469 

(12) q^^ = 0.0422 X |^ = 0.469 
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Member 


Q 


8 


q 


qs 


qs'* 


w lb. /in 


F.S. 


66.7 


11. C4: 




0.475 


5.33 




R.S . 


141.0 


2S.76 


0.0199 








U.S. 


6,01 


4.19 


0.469 


1.965 


8.24 




L.S. 


6.01 


4.19 


0.469 


1.965 


8.24 


0.0594 










Sqe^ = 


33.09 





(13) Wy ^ 0.0594 ^ 0.998 = 0.0595 lb. /in. 

(14) - v/j^ = 0.0572 X 0.0595 = 0.0034 

(15) 11.24 Wp + 23.76 Wj^ = 1 - 4.19 (0.0595 + 0.0594) = 

= 1 - 0.499 = 0.501 

(16) 11.24 Wp + 23.76 (wp - 0.0034) = 0.501 
35 Wp = 0.581;, Wp = 0.0166 lb. /in. 

(17) = 0.0166 - 0.0034 = 0.0132 lb. /in. 

Ip = 2 X (4 X 5^ + 1/12 X 2 X 8) + 1/12 X I/8 x 12^ = 202.6 + 

+ 18 = 220.6 in. * 
Ij^ = 101.3 + 18 = 119.3 in. * 

= II = (4 x 5.57^ + 13.33^ x 2) + 1.33 + 0.17 + l/32 x l/l2 X 

X 21.5^ = 534 + 1.5 + 25.9 = 551.4 in. 

551.4 = 2(1/12 X 21.5^ - 1/12 X A^); 21.5^ - A^^ = 

= 3368 - 9940 - A^ 

A = 18.70 in. = 2 K. 
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4.00 in. 4.00 in. 9.35 in. 9.35 in. 

6.00 in. 6.00 in. 10.75 in. 10.75 in. 

1/4 in. 1/4 in. l/l6 in. l/l6 in. 

2.25 in. 1.25 in. 2.06 in. 2.06 in. 

220.6* in. 119.3* in. 561.4* in. 561.4* in. 

8 X SeVOO^x 220. 6^ x 6^ - 4 ( 6* -| x 6^ x 4^ + 



-8 



+ 1/5 X 4*)] = 1.99 X 10 



^ 2.225^ X 4 ^ 2.25 x 64 ^ 



2 X 86700 X 220.5^ 4 4x3 



2 « -6 



X 20 + ^ X 4 ] = 1.080 X 10 

15 X lb 



3 

¥4i ^ ??q'^5 X a-. = 1.99 X 0.555 x 3.42 X lo" 
o . do 119 . 3 J? 

= 3.77 X 10~^ 



2 X 86700 X 119.3^ 4 '4x3 



X 64 X 20 + -, ^ ^ T ^ X 4^] = 1.243 x 10 
xo X lb 



2«05 f- 8 Y m Q 'z^ fin -^c; * 



8 X 86700 X 561.4 1^ 



[« X 10.75 - 9.35 (10.75*- 



- 2/3 X 10.75 X 9.35^ + l/S X 9.35"^] 
= 4.72 X 10~^ 
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7=7= 1§ ^2.06^ X 9.55 (io,75^ - 9.35^)^ + 

U L 2 X 86700 x 561.4^ 



x 9.35® (10.75^ _ 9.35^) -f 



^ X 9.35^ = 2.62 X 10 ^ 



15 X 256 

= 2/3 (a + 7) w^l"^ = 2/3 x 200^ x (a + 7)w 



= 5,330,000 (a + 7)w® 



W = . w^L^ ^ X 200^ = X 6150 

-D 40 E I 40 X 1.3 X 10^ X I ^ 



I a + 7 Wg Wrp 

F 220.6 1.100 X 10~® 5.87 Wj,^ 27.9 w^^ 33.8 Wj,^ 

R 119.3 1.281 X 10~^ 6.83 Wj^^ 51.5 Wj^^ 58.3 Wp^^ 

U 561.4 2.62 X lo"^ 13.88 11,0 w^^ 24.9 w^^ 

L 531.4 2.62 x lo"^ 13.88 w^^ 11.0 w^^^ 24.9 w^^ 



20 Wj, + 10 = 1 

Wy = 0.10 - 2Wp 

ZW = 92.1 Wj,^ + 49.8 (0.10 - 2 Wp)^ 

= 184.2 w„ - 2 X 2 X 49.8 (0.10 - 2 ?;_,) = 0 
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582.6 Wj, _ 19.9 = 0^ Wj, = 0.0342 

w„ = .0342 lb. /in. 

= .0342 lb. /in. 

= .0315 lb, /in. 

= .0315 lb. /in. 



w. 



-i- = 1 + |3i8 ^ 3 ^ 53.8 = 1 + 0.580 + 2 X 1.357 = 4.294 
Qj, bo.o «3 4.9 

= 0.233 









0.580 X 


0.233 = 0.135 










= 1.357 


X 0.233 = 0.315 






Q 


6 


q 


qs qs ^ 


w 


F 


33.8 


7.34 


0.233 


1.710 12.56 


.0342 


R 


58.3 


12.66 


0.135 


1.710 21.53 


.0342 


U 


24,9 


5.0 


0.316 


1.580 7.90 


.0315 


L 


24.9 


5.0 


0.316 


1.580 7.90 


.0315 










2qs'' = 49.99 






Fig - 2 Sample epar No . 1 
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Cantilever length = 200" 



Scale: 1 inch = 5 inches 



Fig. 3 Sample spar No. 3 



